^A103523 


gONALP^AAVIDSW 


RESEARCH  AND  TECHNOLOGY  DEPARTMENT 


AwroMd  for  public  rcloMt.  dinribution  unlimited. 


NAVAL  SURFACE  WEAPONS  CENTER 

DtMfran.  Virtinis  22448  •  SiNw  Spring.  Mwyland  20910 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  rlFhwi  Data  Enfnd) 


REPORT  DOCUMENTATION  PAGE 


«.  title  rand  Su6U(/e; 

FREE  ELECTRON  LASER  INSTABILITY  FOR  A 
RELATIVISTIC  SOLID  ELECTRON  BEAM  IN  A  HELICAL 
WIGGLER  FIELD 


7.  AUTHORr*; 

Han  S.  Uhm,  Ronald  C.  Davidson,  and 
Science  Applications,  Inc. 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Naval  Surface  Weapons  Center 
Code  R41 

White  Oak,  Silver  Spring,  Maryland  20910 


II.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  RECIPIENT'S  CATALOG  NUMBER 


S.  TYPE  OF  REPORT  «  PERIOD  COVERED 

Final  -  Jan  1981 


6-  PERFORMING  ORG.  REPORT  NUMBER 


S.  CONTRACT  OR  GRANT  NUMBERC*; 


10'  PROGRAM  element.  PROJECT.  TASK 
AREA  A  WORK  UNIT  NUMBERS 


>2.  report  date 


January  1981 


13  NUMBER  OF  PAGES 

43 


U  MONITORING  AGENCY  N  AME  A  AODRESS^I  dlltarani  from  Conirolling  OlllctJ  IS.  SECURI  TY  CL  ASS.  ra<  'Ala  raport; 

UNCLASSIFIED 

ISa.  DECLASSIFICATION  DOWNGRADING 

schedule 


IS.  distribution  STATEMENT  (at  thlt  Raporl; 


Approved  for  public  release,  distribution  unlimited. 


17.  distribution  statement  (of  tho  mbttract  ontorott  In  Block  20.  It  <iHloront  from  Roport) 


19.  KEY  words  (Continuo  on  rmvoroo  oldo  If  nocooomry  mnd  Idontify  by  block  numbmr) 

Free  Electron  Laser 
Microwave  Emission 
Relativistic  Electron  Beam 
Helical  Wiggler  Field 


20.  ABSTRACT  (Continuo  on  rovoroo  oldo  If  nocmoomry  ond  Idontify  by  blocA  numborj 

The  free  electron  laser  instability  for  a  solid  relativistic  electron  beam 
propagating  in  combined  transverse  helical  wiggler  and  uniform  axial  guide 
fields  is  investigated  within  the  framework  of  the  linearized  Vlasov-Maxwel 1 
equations.  It  is  assumed  that  v/vj^  <<1 ,  where  v  is  Budker's  parameter  and 
Yijmc"^  is  the  electron  energy.  Stability  properties  are  investigated  for  the 
choice  of  equalibrium  distribution  function  in  which  all  electrons  have  the 
same  value  of  the  linear  combination  of  transverse  and  helical  invariants, 


DD  I  JAN *73  1473  edition  OF  I  NOV  «5  IS  OBSOLETE 

S/N  0  10  2-014-  «60  I 


UNCLASSIFIED 

security  classification  of  this  page  'IWiAn  OalA  BnlaradJ 


UNCLASSIFIED 


-n-uwiTv  classifiotiom  of  this  p/koermim  d«i«  Entumd) 

20.  -  2YuiniDuCu  =  const.,  and  a  Lorentzian  distribution  in  the  axial 

invariant  c  .  (Here  wu  is  a  constant.)  The  instability  growth  rate  is 
calculated  Tncluding  a°detenni nation  of  the  optimum  value  of  the  ratio  of 
beam  radius  to  conducting  wall  radius  (Ro/Rq)  for  maximum  growth.  It  is 
found  that  tne  maximum  growth  rate  for  a  solid  electron  beam  is  com¬ 
parable  to  that  for  a  hollow  beam  with  similar  parameters.  Moreover,  the 
introduction  of  a  small  axial  momentum  spread  (A/vi^mc  a  few  percent) 
significantly  reduces  the  instability  growth  rate. 


UNCLASSIFIED 


UCUMITV  CLA$IIFICATION  OF  THIS  PAOeOFAwi  OMa  KntararfJ 


NSUC  TR  81-145 


FOREMORO 


^ The  free  electron  laser  instability  for  a  solid  relativistic  electron  beam 
propagating  in  combined  transverse  helical  wiggler  and  uniform  axial  guide 
fields  is  investigated  within  the  framework  of  the  linearized  Vlasov-Maxwell 
equations.  is  assumed  that  v/yjj  «1 ,  where  v  is  Budker's  parameter  and 
Yj,mc  is  the  electron  energy.  'Stability  properties  are  investigated  for  the 
choice  of  equalibrium  distribution  function  in  which  all  electrons  have  the 

same  value  of  the  linear  combination  of  transverse  and  helical  invariants., 

/ 

-  ZYjjtnujjC^  =  const.,  and  a  Lorentiian  distribution  in  the  axial  invariant  C^. 
(Here  is  a  constant.)  ^The  instability  growth  rate  is  calculated  including  a 
determination  of  the  optimum  value  of  the  ratio  of  beam  radius  to  conducting 

CV 
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1.  INTRODUCTION 

1-9 

In  recent  years,  the  free  electron  laser  instability  has  been 
extensively  investigated  with  particular  emphasis  on  applications 
to  Intense  microwave  generation.  For  the  most  part,  previous  theoretical 

analyses  of  this  Instability  have  been  carried  out  either  for  uniform 
3—6 

density  beams  with  infinite  transverse  dimension,  or  for  annular 
7  6 

electron  beams.  ’  The  present  paper  examines  the  Influence  of  finite 
radial  geometry  on  the  free  electron  laser  instability  for  a  solid 
electron  beam  propagating  in  combined  transverse  helical  wiggler  and 
uniform  axial  guide  fields.  The  analysis  is  carried  out  v;lthin  the 
framework  of  the  linearized  Vlasov-Haxwell  eouations.  Including  a 
determination  of  the  optimum  value  of  beam  radius  for  maximum 
growth  rate. 

The  present  analysis  is  carried  out  for  an  infinitely  long  relativ¬ 
istic  electron  beam  propagating  in  the  combined  transverse  wiggler  and 
uniform  axial  guide  fields  described  by  Eq.  (1).  Equilibrium  and  stability 
properties  are  calculated  for  the  specific  choice  of  electron  distribution 
function  (Eq,  (5)], 

-  2Yj^mTi)G(C^)  , 

where  n^,  and  T^  are  positive  constants,  Cj,,  Cj^,  and  are 

the  transverse,  helical,  and  axial  invariants'^  defined  in  Eqs. 

(6)  -  (8),  and  the  axial  distribution  function  is  normalized  according 
to  dC^G(C^)  “  1.  Equilibrium  properties  are  investigated  in  Sec.  II, 

and  stability  properties  are  examined  in  Secs.  Ill  and  IV,  assuming 

2 

that  v/yjj  «  1,  where  v  is  Budker's  parameter,  Yi^mc  is  the  characteristic 

electron  energy,  and  m  is  the  electron  rest  mass.  j 
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In  Sec.  Ill,  making  use  of  the  linearised  Vlasov-Maxwell  equations, 
we  obtain  the  coupled  eigenvalue  equations  (33),  (37),  (40),  (45),  and 
(46)  that  describe  free  electron  laser  stability  properties  In 
circumstances  where  the  perturbed  transverse  fields  can  be  approximated 
by  the  vacuum  waveguide  fields.  For  short  wavelength  perturbations, 
the  axial  component  of  the  perturbed  longitudinal  field  can  be  approximated 
by  [Eq.  (49)], 


\  0  ,  ot 


r/Ro)  . 


0  <  T  < 


otherwise 


where  Jj^(x)  is  the  Bessel  function  of  the  first  kind  of  order  t, 

6,  ,  is  the  s'th  zero  of  J,(6,  ,)  -  0  and  <>  ,  is  a  constant.  In 

Sec.  IV,  substituting  Eq.  (49)  into  the  coupled  eigenvalue  equations, 
we  obtain  closed  algebraic  dlspersi02i  relations  foj  the  transverse 
electric  (TE)  and  transverse  magnetic  (TM)  polarizations. 


Introducing  the  normalized  dimensionless  function  [Eq.  (59)], 

2 

‘  ^ei,s-  7,27  J-72  ’ 

it  is  shown  in  Sec.  IV  that  the  coupling  between  the  longitudinal  and 


transverse  perturbations  Is  proportional  to  G  t  (a  . ,  ^ 

£ S  Et  X  g  8  U  C 

for  the  TE  mode,  and  to  G.  i(6.r<i  )  for  the  TH  mode.  Here 

is  i+l,s  0  c 

“i+l.s  “***  ^+l,s  sth  zeroes  of  -  0  and  “  0, 

respectively,  R^  is  the  radius  of  the  outer  conducting  wall,  and  the  prime 

denotes  (d/dx)J  (x).  Assuming  that  the  maximum  of  the  function  C  , (x) 
i  £s 

occurs  at  X  ~  X  , ,  we  note  that  the  maximum  instability  growth  rate 
is 

occurs  at  a  value  of  R_/R  given  by  R^/R  -  x.  for  the  TE 

0  c  0  c  is  i+l,s 

mode,  and  by  Rq/R^  ”  s  ^or  the  TM  mode.  This  result  is 

g 

different  from  that  obtained  for  a  hollow  electron  beam. 
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A  detailed  numerical  analysis  of  the  TE  mode  [Eq.  (62)]  and  TM 

node  [Eq.  (63)]  dispersion  relations  is  presented  in  Sec.  IV.  Two 

features  are  noteworthy  from  the  numerical  analysis.  First,  for  the  optimized 

value  of  Rn/R  ,  the  instability  grovith  rate  for  the  TH  mode  is  comparable 
u  c 

to  that  for  the  TE  mode.  Moreover,  the  growth  rate  is  reduced  substantially 
by  introducing  a  small  amount  of  axial  momentum  spread  (A/y^nc  0.01). 


9/10 
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II.  EQUILIBRIUM  THEORY  ANU  BASIC  ASSUMPTIONS 


The  equlllbrluni  configuration  consists  of  a  relativistic  electron 
beam  propagating  in  the  combined  transverse  wiggler  and  uniform  axial 
guide  fields  de.scribed  by 


B  =  -«Bcos(6  -  6Bsln(e  -  k^zle^  B^e^ 


where  Bq  and  6B  are  constants,  and  kg  is  the  axial  wavenumber  of  the 
helical  wiggler  field.  In  Eq.  (1),  cylindrical  polar  coordinates  (r,e,z) 
are  used,  with  z-axls  along  the  propagation  direction,  and  and 

are  unit  vectors  in  the  r-,  6-,  and  z-directions ,  respectively.  In 
the  present  analysis,  we  assume  that  tlie  axial  wavenumber  of  the  helical 
wiggler  field  is  sufficiently  large  that 

T~---  I-  “I  «  UpU  ,  (2) 

'“O  ■  “c'  ^0^  ° 

where  ”  *^o'^b’  ^^c  electron  cyclotron  frequency, 

2 

Rq  is  the  characteristic  beam  radius,  is  the  characteristic 

2  1/2 

electron  energy,  c  is  the  speed  of  light  in  vacuo,  ^’^b 

is  the  mean  axial  velocity  of  the  electron  beam,  and  -e  and  m  are  the 
electron  charge  and  rest  mass,  respectively. 

It  is  also  assumed  that 


v/yk  « 


where  v  “  N.e  /me  is  Budker's  parameter. 


,2x  rR. 

■  J„  *  1 


dr  r  nj^(r,  6  -  k^z)  , 


is  the  number  of  electrons  per  unit  axial  length  of  the  beam,  nj^(r,  6-  kjjZ) 
is  t'tie  equilibrium  electron  density,  and  R^.  is  the  radius  of  the  conducting 
wall.  The  inequality  in  Eq.  (3)  indicates  that  the  beam  is  very  tenuous, 


11 
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In  the  parameter  regimes  of  practical  interest  for  free  electron 

laser  applications,  the  axial  distribution  function  G(C^)  is  strongly 

peaked  about  »  ^b'"'^b  "  const.  ,  with  characteristic  half-width 

LC  «  Y.mV.  .  Moreover,  in  the  present  analysis.  It  is  also  assumed 
Z  D  D 

that  the  axial  motion  is  nonresonant  with 


1.2  2  ,  2 
*  “c  • 


where  axial  velocity  of  a  typical  beam  electron. 


(10) 


We  therefore  approximate  Eq.  (8)  by 


iti  iE 


10 


0  c  0 


(11) 


Making  use  of  Eqs.  (6),  (7),  and  (11),  it  is  straightforward 


to  show  cliat  the  combin.U  i<>o  Cj,  •  In  Kq.  (5)  i-hu  be  expressed  as 


iO 


Cx  -  2rj^mm^Cf,  -  (p^  -  )“-(e  -  k(j*)J 

d 

+  (pQ  -  Yb^-V  ^  -  M>) 


(12) 


+  Yj<n%o<’'*  ®  ’ 


where  idq  -  *^0^’  “c  '  effective  potential  il'pir,  6  -  k^z) 

is  defined  by 


/"“c  ~  r  6E 
"c'%  ^“0  ■  “c/  *^0  ®0 

■  ( !:g-'J!by  <  / M  f2  ^ 

\“0  ■  “c/  Wt  ®o'  V\  “0  /  ^  ^ 


i2  2 


sin(e  -  kgZ) 


(13) 


As  a  simple  example,  we  consider  an  axial  distribution  function 


in  which  all  electrons  have  a  same  value  of  C_,  i.e., 

z 
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where  the  mean  radius  is  defined  by 


y,  m  \  U-  -  u  /  2  \  B.  /  b 

h  0  c  0  J 


<“b“c  ■  “b^ 

and  use  has  been  made  of  Eq.  (2).  Additional  general  equilibrium 
properties  associated  with  tiie  ilistribution  function  in  Eq.  (3), 


including  helical  distortions  of  the  beam  equilibrium  for  finite  6B/Bq, 
are  discussed  in  Ref.  10. 


14 
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111.  LINEARIZED  VLASOV-MAXWELL  EOLATIONS  FOR  A  TENUOUS  BEAM 


In  this  section,  we  make  use  of  the  linearized  Vlasov-Kaxwell 
equations  to  investigate  the  free  electron  laser  stability  properties 
of  a  relativistic  solid  electron  beam  described  by  the  equilibrium 
distribution  function  in  Eq.  (5).  We  adopt  a  normal-mode  approach  in 
which  all  perturbations  are  assumed  to  vary  with  time  and  space 
according  to 

6i(x,t)  =  I  i*i^"\r)exp{ll£6  -E  (k  -t  nkjj)z  -  ut]}  ,  (17) 

i,n 

where  Imui  >  0.  Here,  u)  is  the  complex  eigenf  requency  ^  k -E  nkp  is  the 
axial  wavenumber,  and  t  and  n  are  integers.  Moreover,  it  is  also  assumed 
that  the  perturbations  are  close  to  resonance  with 


(uj  •  (k  +  nkQ)V|^|  ,  (18) 

where  =  k^V^^  and  ui^  «  eBQ/yj^mc. 

The  Maxwell  equations  for  the  perturbed  electric  and  magnetic  field 
amplitudes  can  be  expressed  as 


^  X  i(x)  »  i(m/c)g(x)  , 

V  X  B(x)  -  (4ii/c)i(;|i)  -  i(u./c)E(x)  , 


(19) 


where 


is  the  perturbed  current  density. 
fj^(x,ig)  »  e  di  exp(-iut) 


In  Eq.  (20), 


f<?’>  + 
%  % 


v'  X  B(x') 


(20) 


(21) 
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is  the  perturbed  distribution  function,  i  -  t*  -  t,  and  the  particle 
trajectories  ;(c'(t')  and  ]j'(t')  satisfy  d;|c'/dt'  -  and  d;g'/dt'  “  -e;^'  *  ^/c, 
with  "initial"  conditions  •  t)  '  and  v'(t*  «  t)  • 

Within  the  context  of  Eqs.  (3)  and  (IB),  the  perturbed  distribution 
function  in  Eq.  (21)  can  be  approximated  by 


fj,(;^,jg)  ~  ~  ~  j  exp(-ian)|2^Ymiui(v'  •  ]^) 

^  •  h  ?)!?7  *b}  ■ 


(22) 


2  22  1/2 

where  y  “  (1  +  jj  /m  c  )  ,  and  use  has  been  made  of  Eq.  (19).  To  lowest 


order,  the  axial  motion  of  an  electron  is  free-streaming  with 


10 


ym 


t)  . 


(23) 


Moreover,  within  the  context  of  Eq.  (18),  on  the  right-hand  side  of 
Eq.  (22)  we  retain  contributions  to  and  in  the  orbit  integral  of 
the  form® 


6  B 

«  y  - -  cos(e  -  k^z  -  k^v^r)  , 

r  Z(4)«-wB«  0  Oz  ’ 

0  c  U 


(24) 


and 


^  "O  "  “c  ®0 


^  sin(e  -  kQZ  - 


(25) 


Finally,  since  the  oscillatory  modulation  of  the  radial  and  azimuthal 
orbits  is  small  amplitude  [Eq,  (2)],  we  approximate 


I.-'  =  r  ,  e’  =  6  , 


(26) 


in  the  arguments  of  the  perturbation  amplitudes  on  the  right-hand  side  of 
Eq.  (22). 
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Substituting  Eqs.  (23)  -  (26)  Into  Eq.  (22),  we  obtain  the  perturbed 
distribution  function 


^  f^"^exp{l[i0  -1  (k  -1  nkQ)z)) 


l.n 


lec  I  (k^nkp)zl)/ 

“  t,n  “  -  +  ">‘0>''z  \  “  * 


(r)  + 


e6B 


2Yi!ic^kj|  ***0 


(27) 

Wiiere  the  function  (jg.iu.k)  is  defined  by 


afj  3f? 

An'^lf*^**^)  '  2(Yinu  -  (k  -1  ”'kQ)p^]  — 2  +  (k  +  "'•‘q)  ^ - 


(28) 


3  Pi 


and  g  •  V  /c.  In  Eq.  (27),  the  term  proportional  to  X  is  the 
z  2  n 


longitudinal  portion  of  the  perturbed  distribution  function. 

Similarly,  the  terms  proportional  to  X^^j^  and  X^_j^  in  Eq.  (27)  are  the 
transverse  electromagnetic  portions  of  the  perturbed  distribution  function. 


Consistent  with  Eq.  (18),  the  eigenf requency  ui  can  be  approximated 


2  2  2 

by  u)  =  (k  +  nkQ)Vj^.  We  therefore  approximate  u  /c  -  (k  +  nkp  -f  k^)  by 


u^/c^  -  (k  -f  nkg  +  * 


(k  +  nk.) 

- +  2kjj(k  +  nkj,)  -t  k^ 


<  0  , 
(29) 


for  k  +  nkg  >  0.  Evidently,  Eq.  (29)  indicates  that  the  n  -t  1  mode  in 
Eq.  (27)  is  a  non-propagating  wave  in  a  vacuum  waveguide.  Without  loss 
of  generality,  for  a  tenuous  beam,  we  therefore  assume 


p(n-H) 

r.l-1 


(r) 


f(n-H) 

6,1-1 


(r) 


0  . 


(30) 


in  the  subsequent  analysis.  Making  use  of  Eq.  (30),  f"  Eq-  (27) 

can  then  be  expressed  as 


17 


NSWC  TR  81-145 


jgj,  exp{l(ie  +  (k  +  nkplz)) 

“  u  I  u  -  (k  +  nk.)v 

£  ,n  u  z 


X  B  E^"J(r) 

I  D  Z  Z  f  X> 


!• 


(31) 


(32) 


where  the  dloenslonless  parameter  A  is  defined  by 

eAB  *^0 

,  2,  u-  -  a.  * 

2Yj.mc  kg  0  c 

and  use  has  been  made  of  the  approximation  Y  ~  which  is  consistent 
with  Eq.  (18). 

From  Poisson's  equation,  7  •  j|(ijj)  “  ^rp(x).  and  the 
Maxwell  equation  (19),  we  obtain  the  differential  equation, 


[^,2  ^  .  (k  +  nk^)2]i<';i 


4xi(k-^nko) 

(r)  =  - ^ -  Pj^  (r) 


(33) 


for  tlie  axial  ( I  jn^iitudinal)  component  of  the  perturbed  electric 


field  E^"^.  In  Eq.  (33),  p^^^r) 


-e  ( 


d^pf^”^  is  the  perturbed  charge 


density,  7?  i  r"^0/3r)  (r3/3  r)  -  £^/r^,  and  use  has  been  made  of  J^"^(r) 
V|^pJ^"^(r).  In  the  tenuous  bean  limit  lEq.  (3)],  the  transverse  field 
components  in  Eq.  (31)  can  be  approximated  by  the  vacuum 

O 

waveguide  fields.  In  this  context,  the  present  stability  analysis 
utilizes  the  vacuum  transverse  electric  (TE)  and  transverse  magnetic 
(IM)  waveguide  modes  as  a  convenient  basis  to  represent  the  general 
electromagnetic  field  perturbation  which  is  determined  from^^ 


^  -  (k  +  nkjj 


ko) 


(3«) 


’  viz  ^z"t+i^^^  ^  c  ?z  *  ^^®z,£+l^^^ 


i(n-l), 


Making  use  of  Eqs.  (19)  and  (34),  and  neglecting  the  perturbed 
current  density,  the  vacuum  waveguide  fields  can  be  expressed  as 
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for  the  TE  mode,  and 
r(n-l) 


(oR 


(35) 


^“.4.,  r/R  )  . 


“t+l.s 


2,1+1 

;.(n-l) 


■  ^l+l.sVl<^-H,s^/«c>  • 


(36) 


F  — f.x  _  '”‘0  -  V«c 

r.l+1^  i  -  4+l.s"l(ei+l.s’^/«c>  • 

for  the  m  »ode.  In  Eqs.  (35)  and  (36).  and  are  constants. 

Jj^.Cx)  is  the  Bessel  function  of  first  kind  of  order  I*  and  a 

i+l,s 

°f  ■';+iW+i,s>  =  0  =  “• 

respectively.  Here  the  prime  (')  denotes  Jj,,(x)  =  (d/dx) (x) . 

Aftei  ouuc  straightfoiwaio  algebraiv,  manipulation  of  Eqs.  (19),  (35),  and 
(36),  it  can  be  shown  that 


f  2  C.2 

l!i tv  J.  _  ,.  \2  1+1, s 

1^2  -~X' 

c 


b  J  f  ^  \ 
l+l.s'^l+ll  R  / 


(37) 


for  the  TE  mode,  and 


^  4  -  (k  +  nk.  -  k  ,)^  -  r  J  (  ) 

c^  0  0  1,2  ‘^l+l.s'^l+lV  R  } 


I 


(38) 


“  ^  "'^0  -  -  7  . 

for  the  TM  mode.  Moreover,  making  use  of  the  continuity  equation, 

0  •'O^-'z.l+l  r  ar  '’"‘^r.l+lJ  ^  r  '^6,1+1’ 

(39) 
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the  approximation  (consistent  with  Eq .  (3)], 
and  approximating  k  -f  nk^  =  k^/ (1  -  V^/c)  on  the  right-hand  side  of 
Eq.  (38),  we  find  that  Lq.  (38)  can  be  expressed  as 


^  -  (k  +  nkjj  - 

c 


ko) 


e+l ,  s 


E  J  f!i±ir£ 

t+l.s'^t+U  R 


An  Is _ 

rc  |3r 


[rJ 


(n-1) 

r,£+l 


(r)]  +  i(t+l)3J’;;ii(r)}  . 


(40) 


for  the  TM  mode. 

For  convenience  of  notation  in  the  subsequent  analysis,  we  intro¬ 
duce  the  effective  susceptibility. 


Xn  n>  d^p 


»  tiJ  *  ^) 

w  *  (k  +  nkQ)v^ 


(41) 


^'oreover,  to  simplify  the  present  analysis,  we  also  assume  tliat  the 
beam  rotation  is  slow  with 


“b  “c’  “O  • 

Within  the  context  of  Eq.  (42),  we  can  show  from  Eq.  (12)  that  the 
equilibrium  distribution  function  is  an  even  function  of 


Pj  -  2Y^mcAcos(e  -  kpZ)  (43) 

and 

Pg  +  2YjjmcAsin(e  -  k^z)  ,  (44) 


for  the  beam  rotations  satisfying  <<  Wq.  Making  use  of 
Eqs.  (31),  (41),  and  (43),  the  perturbed  charge  and  current  densities 
are  given  by 


;(n-l) 

e,i+i 


(r) 


Anu 


Gj^(w,k,r) 


- (n-1)  ,  . 
r,t+l'  •’ 


icAp^"^ (r) 


(45) 


where  the  function  Gjj(u,k,r)  is  defined  by 


20 


NSWC  TR  81-145 


CjjCw.k.r) 


»  6 
,  n  z  z  ,i. 


(r) 


-(n-1)  .  .  . 


(A6) 


Equations  (33),  (37),  and  (40),  when  combined  with  Eq.  (45),  constitute 
one  of  the  principal  results  of  this  paper  and  can  be  used  to  investigate 
stability  properties  for  a  broad  range  of  system  parameters.  Moreover, 
in  limiting  cases,  the  dispersion  relation  for  the  free  electron  laser 
instability  can  be  obtained  in  a  closed  form  (Sec.  IV). 
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IV.  FREE  ELECTRON  LASER  STAEILITY  PROPERTIES 


In  this  section,  simplified  expressions  are  obtained  for  the  longi¬ 
tudinal  perturbations  in  Eq.  (33), and  the  results  are  used  to  derive  the 
dispersion  relation  for  several  values  of  azimuthal  harmonic  number  t. 

The  present  analysis  assumes  short  wavelength  perturbations  with 

|q^|  =  |(k  -*■  nkp)^  -  Jlc^\  »  1/r2  .  (A7) 

Moreover,  for  <u  =  (k  -f  nkQ)V^  and  k  +  nV.^  =  kp/{l  -  V^^/c),  the  inequality 
in  Eq.  (47)  can  be  expressed  in  the  equivalent  form, 

(1  ^ 


which  is  readily  satisfied  in  the  parameter  regimes  of  present  experimental 
interest.  As  shown  in  Appendix  .\,  for  short  wavelength  perturbations 
satisfying  Eq.  (48),  the  axial  component  of  the  perturbed  electric  field 
E^''j(r)  in  Eq.  (33)  can  be  approximated  by 

Z ,  I 


c(n) 


(r) 


{^l,s''^E^^£.,s 


r/R^j)  ,  0  <  r 

otherwise  . 


<  R„ 


(49) 


In  Eq.  (49),  6,  ,  is  the  s’th  root  of  J«(C,  ,)  =  0,  and  ,  is  a 

constant. 

Substituting  Eqs.  (35)  and  (49)  into  Eq .  (46),  multiplying  Eqs. 

(33)  and  (37)  by  '^'^l+l^“n-l 

and  integrating  from  r  •  0  to  r  »  R^,  we  obtain  two  homogeneous  equations 
relating  the  perturbation  amplitudes  ,  and  b  .  For  the  TE 

mode  polarization,  these  are 

‘...■t'  '  *<'0  -  '>  (7  ) 

^  °  (50) 
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and 


(k  +  nk.)R 
_  u  c 

'*'b'’£+l,s 


£+1 


fR 

c 

•"Jo 


dr  r 


^  -  (k  +  nk„  -  Vy  -  -h 


'n,n-l 


,2  l°l+l,s^ 


£-H' 


c 

(51) 


where  •(x)  is  the  Heaviside  step  function  defined  by 


«(x) 


f  1  .  X  >  0  . 

[  0  ,  otherwise 


(52) 


Similarly,  for  the  Ttl  mode  polarization,  we  obtain 

,2 


and 


fR 

••lo‘ 


(53) 


»2  /  .  » &  \ 

£+l\  R  ) 
c 

(54) 

s*  '  t  s'^  \ 

(k  +  nkQ)RQ  "Jr  r  «(Ro  '  '^^’‘n.n'^i-H  {  R^  ) 


^  -  (k  +  nk^  -  kj,)"  - 


“  +  A  X, 


n,n-l 


’‘\+l(  P. 


^£-»-l,s^ 


)• 


where  use  has  been  made  of  u  =  (k  +  nk-)V.  and  (k  +  nk  )  =  k,./(l  -  V, /c). 

Ob  Q  U  b 

In  the  present  analysis.  It  is  assumed  that  the  axial  distribution 
function  has  the  form 


G(C  )  =  ^ 


(C,  -  yb"’''b)  A 


2  ’ 


(55) 


where  A  is  the  characteristic  spread  in  C  about  the  mean  value  C  >  r,  mV,  ■ 

z  z  b  b 
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We  further  assume  that  the  characteristic  spread  fi  is  small  in  comparison 
with  Substituting  Eqs.  (5)  and  (55)  into  Eqs.  (28)  and  (41), 

we  obtain  the  approximate  expression 


-  (k  nkQ)(k  -*■  n'kQ)c^ 


n,n 


0  ,  Sq  <  r  _< 


0  <  r  c  R^. 

(56) 


In  obtaining  Eq.  (56),  use  has  been  made  of  Lo.  (18).  Making  use  of  the 

definition  of  Budker's  parameter  in  Eqs.  (3)  and  (4),  the  term  4v/y,R^ 

b  0 

*>222 

in  Eq.  (56)  can  also  be  expressed  as  4v/y,  R)!  =  lo  , /c  ,  where  u  , 

b  U  pb  pb 

2 

4iinQe  is  the  plasma  frequency-squared. 

Tlie  condition  for  a  nontrivial  solution  to  Eqs.  (50)  and  (51) 

is  that  the  determinant  of  the  coefficients  a.  ,  and  b  ,  be  equal 

t,s  t+l,s 

to  ^ero.  After  some  algebraic  manipulation,  we  find  that  the  TF.  mode 
dispersion  relation  can  be  expressed  as 


u  -  (k  +  nkQ)Vj^  +  i 


|k  4  nk^lA 


\4- 


j(  |k  4  nk„|A 

^  .2  \ 

(  “  -  (k  4  nkQ)Vj^  4  i  3  0 

-  4  S 

3„2  / 

U  ^b“ 

^b’^O  J 

4A 


A  ^ 

2  VC 


0  c 


kjj(k  4  nkg  -  kjj) 


•t41,s 


1 


P. 


f °841,s”0 
■‘iss'V 


^"tss-(^T^)h-  (•^  +  nko)V^4i 


|k  4  nkpli 


3 

Yb” 


3„2  iss'V 


*^£41,8^0 


where  the  coupling  coefficient  i (a, . ,  Rn/K  )  is  defined  by 

£SS  y S  V  C 


«tss'(*> 


*t4l.s 


*£41, s 


(i  ^  1)' 


and  the  functions  G^^,(x)  and  defined  by 


(57) 


(58) 
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2  Ji(x) 

ri— :2 


and 


(x" 

,2 


j;^j(x)  -  j^(x)j^^^(x) 

G^^.(x) 


(59) 


(60) 


In  Eq.  (57),  the  subscript  s  and  s'  represent  the  radial  node  numbers 
of  the  transverse  and  longitudinal  perturbations,  respectively. 

For  small  wlggler  amplitude  (A  «  1),  vc  Investigate  free  electron 

laser  stability  properties  for  u  and  k  +  nkQ  near  the  simultaneous 

2  2  2 

zeros  of  the  transverse  dispersion  relation,  u,  -  (k  -f  nk^  -  k^)  c  - 

2  2  2 

°J.-H  s"^  ^^'c  "  longitudinal  dispersion  relation 


u  -  (k  +  nkQ)Vj^  -i  i 


|k  +  nkjjlA 


Y.m 


-  4 


vV 

^b’^0 


0  . 


(61) 


In  this  regard,  making  use  of  Eq.  (61),  the  TE  mode  dispersion  relation 
in  Lq.  (57)  can  be  approximated  by 


i  HL-  -  (k  4-  nk  -  k  )^  -  — -4—  ^ 

12  '  0  ^0''  2  I 

I  c  R  J 

c 


k  +  nk 


u  -  (k  +  nkQ)Vj^  +  i 


i!^r 
j 


-  A 


3„2 


4ft 


2  vc‘ 


b  c 


ko(k  nkjj  -  kjj)  2 


E  /  °t-t-l,s‘^0 
^tss'V  R 


In  a  similar  manner,  from  Eqs.  (S3)  and  (54),  ve  obtain  the 


c 

(62) 


approximate  TIi  mode  dispersion  relation, 

„2  \ 


2 

2  ■ 

-  4 


6 


(k  +  nkp  -  kjj) 


i-*-l.s\ 


-  (k  +  nkQ)V|^  +  i 


|k  +  nkplft 


-12 


2  •) 


3=2 


4ft 


2  VC 


kQ(k  +  nkjj  -  kjj) 


^£-H,s 


I 

^tss’  ( 


(63) 

^£+l,s*^0 


where  the  TK  mode  coupling  coefficient  Q,  , (6. . ,  Rn/®  )  defined  by 

”  £ss  £+l,s  U  c 

(64) 


<‘ss'<’‘>  '  «£s'<*>/42^^41.s>  • 


and  the  function  G^^,(x)  is  defined  in  Eq.  (59). 
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Figure  1  shows  plots  of  versus  x  obtained  from  Eq.  (59) 

for  (a)  s'  “  ®o  1  s'  '  3'  case 

B.  ,  =  Bn  I )  the  plots  of  G  , (x)  for  arbitrary  6„  i  are  similar 

to  those  for  6^^  s'  "  3  1(b).  As  shown  In  Fig.  1(b),  the 

quantities  C*  ,  =  G.  i (x.  ,)  and  x.  ,  denote  the  maxiraum  value  of 
is  is  is  is 

Gj^^'(x)  and  the  corresponding  value  of  x  for  a  specified  6^^  .  For 

example,  in  Fig.  1,  (Xj^^,,  G*^,)  =  (0,  0.69)  for  s'  "  ®0  1 

(Xo.i  .  gJ  , )  =  (9.8,  0.06A)  for  E  ,  «  6,  ».  Shown  in  Fig.  2  are 
tS  *  s  -I  *  J 

plots  of  (a)  ’  corresponding  values  of  G^^,  = 

^ts'  ^  several  different  values  of  the  azimuthal  and 

radial  mode  numbers  i  and  s'.  It  is  evident  from  Fig.  2(b)  that 
gJ^,  decreases  rapidly  with  increasing  values  of  the  mode  numbers  t 
and  s'.  Moreover,  we  note  from  Fig.  2(a)  that  Xj^i  can  be  approximated  b 


’‘is-  ^  Sf.s'  •  ^ 


(65) 


In  this  regard,  for  s'  y  1,  cj  ,  can  be  approximated  by 

US 

^ts'  -  °is'(6i,s'^  •  "ll(^,s'>'  ^ 

r 

Shown  in  Fig.  3  are  plots  of  (a)  Q  , /G  ,  for  the  TE  mode 

iss  £S 

and  (b)  Q£gs''^^£s'  ^  mode,  obtained  from  Eqs.  (58)  and  (64) 

respectively.  Note  that  the  curves  in  Fig.  3  are  independent  of  the 

longitudinal  radial  mode  number  s'.  Evidently,  tlie  values  of  Q  ,/G  , 

iss  is 

and  Q  ,/G  ,  Increase  with  Increasing  values  of  azimuthal  and 

iss  is 

transverse  radial  mode  numbers,  i  and  s.  After  careful  examination 
of  Eqs.  (58)  and  (64),  we  find  that  the  maximum  coupling  between  the 


transverse  and  longitudinal  modes  occurs  for  a  value  of  Rg/Fj,  given  by 
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ts' 

1 

i+l,s 

,  /e 

ts' 

'  '^t+l,s 

V"c  ■  ^ 

Equation  (67)  Is  valid  only  when  ^  for  the  TE  mode,  and 

*is>  ™  ^s-  "  “m.s  *is'  " 

(TM),  the  naxlmun  coupling  occurs  for  Rq/R^.  “  !•  The  maximum  coupling 

coefficients  corresponding  to  Eq.  (67)  can  be  determined  from  Figs.  2(b) 

and  3.  For  example,  for  (t,s,s')  -  (3,2,1),  we  determine  that  the 

maximum  coupling  coefficient  and  the  corresponding  ratio  R^/R^,  are 

given  by  (Q^ggi  Rq/^^c^  “  (1-607,  0.625)  for  the  TE  mode,  and  (Q^ggii  Rq/Rc) 

(1.83,  0.52)  for  the  TM  mode. 

It  is  instructive  to  examine  the  present  results  for  perturbations 

with  the  lowest  mode  numbers,  i.e.,  (I!.,s,s')  »  (0,1,1),  particularly 

for  a  beam-filled  waveguide  with  Rq/Rc.  ”  1-  !•'  this  limit,  from  Fig.  1(a), 

we  obtain  G,  ,  (a  )  »  0.4  for  the  TE  mode,  and  G,  ,(iS,  ,)  “  0.045  mode, 
ts  1,1  Is  1,1 

We  therefore  conclude  that  the  TE  mode  polarization  is  the  most  unstable. 
Multiplying  G  ,  (u,  ,)  -  0.4  by  Q^,,_t/G  ,  -  4.2  in  Fig.  3(a),  the 

£S  X « A  KSS  K.S 

E 

coupling  coefficient  is  given  by  “  1-7.  Assuming  zero  axial 

2  2 

momentum  spread  (a  0)  and  short  axial  wavelengths  (RgR^  D* 
the  TE  mode  dispersion  relation  in  Eq.  (62)  can  be  approximated  by 


2 

(  -  (k  +  nk^  -  k^)  1 

<(u  -  (k  +  nkp)!?!^)^ - 3 

'•  c  ' 

1  ^b  J 

,2  2  k2 

'  “pb^O  ’ 


(68) 


for  the  (1,8, s')  =  (0,1,1)  perturbation  and  R_/R 

0  c 


1.  Equation  (68) 


is  similar  in  form  to  the  result  obtained  by  Davidson  and  Uhm  for 


a  uniform  density  beam  with  Infinite  cross  section.  In  particular, 
the  constant  numerical  factor  on  the  right-hand  side  of  Eq.  (68)  is 
equal  to  3.4,  whereas  in  Ref.  3  the  constant  numerical  factor  is  equal  to  8. 
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Finally,  we  have  investigated  detailed  stability  properties 
by  solving  the  dispersion  relations  in  Eqs.  (62)  and  (63)  numerically 
for  a  broad  range  of  system  parameters.  Defining  the  normalized 
Doppler-shifted  eigenf requency  by 

n  -  (<u  -  (k  -f  nkQ)Vj^)/kpC  ,  (69) 

we  calculate  the  normalized  growth  rate  =  lm!i  from  Eqs.  (62)  and 

(63).  Shown  in  Fig.  4  are  plots  of  the  normalized  growth  rate 

versus  (k  -F  nk^l/k^  for  (t,s,s')  =  (3,2,1),  k^^R^  =  10,  '  10,  = 

2 

0.02,  and  A  “  0.01,  with  (a)  R-/P.  =  x,,/a,  .  for  the  TE  mode,  and 

U  c  Ji  4,2 

(b)  R„/R  =  x.,/6,  for  the  TM  mode.  For  these  optimized  choices  of 
U  c  31  4,2 

Rfl/Rci  the  instability  growth  rate  for  the  TI-I  mode  is  comparable  to 

that  for  the  TF.  node.  Moreover,  the  growth  rate  is  reduced  substantially 

by  introducing  a  small  amount  of  axial  momentum  spread  (i/y^^mc  0.01). 

We  conclude  this  section  by  pointing  out  two  areas  in  which  the 
analysis  can  be  extended.  First,  the  restriction  to  very  short  |Vave- 
length  perturbations  [Eq.  (46)]  can  be  removed  in  a  relatively  straight¬ 
forward  manner.  Second,  p.arallelliig  the  self-consistent  theoretical 

g 

formalism  developed  in  previous  studies,  the  stability  analysis  can 
also  be  carried  out  without  making  the  approximation  that  the  transverse 
perturbations  are  represented  by  the  vacuum  waveguide  fields. 


29/30 


NSWCTR  81-145 


V.  COSCLLSIONS 

In  this  paper,  we  have  examined  the  free  electron  laser  instability 

for  a  solid  relativistic  electron  bc2un  propagating  in  the  combined 

transverse  uiggler  and  unifonr.  axial  guide  field  given  in  Eq.  (1). 

The  analysis  V7as  carried  out  within  the  framework  of  the  linearized 

Vlasov-Maxwell  equations.  The  equilibrium  (Sec.  II)  and  stability  (Secs. 

Ill  and  IV)  properties  were  Investigated  in  detail  for  the  choice  of 

distribution  function  in  wiiich  all  electrons  have  the  same  value  of 

the  linear  combination  of  transverse  and  helical  invariants,  C,  -  2y,niuj,  C,  , 

b  b  h 

and  a  Lorentzian  distribution  in  the  axial  invariant  [Cqs.  (5)  and 

(18)].  One  of  the  most  important  conclusions  of  this  analysis  is  that 

the  maximum  instability  growth  rate  for  a  solid  electron  beam  is  comparable 

to  tliat  of  I  iioilow  be.nn  vifi  siifilai  p.'irami.'tei  s.  Kureover,  it  is  also 

found  that  the  maximum  growth  rate  occurs  at  a  value  of  corresponding 

to  Rn/R,.  •  t/a.,,  for  TE  mode  perturbations,  and  R,,/R  «  x, 

^  c  £s'  t-H,s  ’  0  c  is  f-H,s 

for  TM  mode  perturbations.  For  these  optimized  values  of  R,.,/R  ,  the 

U  c 

instability  growth  rate  for  the  Tfi  mode  is  comparable  to  that  for  tlie  TE 
mode.  Moreover,  the  growth  rate  is  substantially  reduced  by  introducing 
a  small  amount  of  axial  momentum  spread  (i/yj^mc  0.01). 
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FIGURE  2  PLOTS  OF  (a)  Xg3'  AND  (b)  THE  CORRESPONDING  G^^’  =  Gg^'  (xg^')  FOR  SEVERAL  VALUES 
OF  AZIMUTHAL  AND  RAIDAL  MODE  NUMBERS,  £  AND  S'. 
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(b)  TM  mode 
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FIGURE  3  PLOTS  OF  (a|  Q^ss'/Ggs'  <•»> 

EOS.  (58)  AND  (64)  FOR  SEVERAL  VALUES  OF  K  AND  S. 
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FIGURE  4(a)  PLOTS  OF  NORMALIZED  GRWTH  RATE  Slj  VERSUS  (k  +  nko)/ko  FOR  IB.S.s')  -  13,^1),  7^  “  ^0, 
Vh^  -  aoz,  AND  »  aoi,  with  (a)  Rq/R^  -  *31/04,2  POR  the  te  mode. 


FIGURE  4<b)  Rq/Rc  =  2  THE  TM  MODE 
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APPENDIX  A 

LONGITUDINAL  PERTURBATIONS  FOR  THE  FREE  ELECTRON  LASER  IH'^tartLITY 


In  this  Appendix,  we  investigate  properties  of  the  longitudinal 
perturbations  about  an  electron  beam  propagating  through  a  cylindrical 
waveguide  with  radius  R^.  In  the  present  analysis.  It  Is  assumed  that 
the  perturbations  have  short  wavelength  with 

|q^|  -  |(k  -E  nkj,)^  -  »  1/Rq  ,  (A.l) 

which  can  also  be  expressed  as 


for  the  frequencies  of  interest  for  free  electron  laser  applications. 
Equation  (A. 2)  is  easily  satisfied  in  parameter  regimes  of  present 
experimental  Interest.  In  the  limit  of  a  small  wlggler  amplitude  (A  -*  0) , 
we  obtain  the  longitudinal  eigenvalue  equation, 


/  1  3  i_  2  \  ;(n)  ,  , 

(  r  3r  3r  ■  2  "  %  ) 


3, 2  ’ 


(w  -  (k  -f  nkglVj^  +  i(k  +  nkplfi/r^m] 
from  Eqs.  (33),  (A5),  (46),  and  (36).  In  Eq .  (A. 3),  €(x) 


step  function  defined  in  Eq.  (52),  and  u 


pb 


Avc^/^bl^S  is 


frequency-squared. 

For  notatlonal  simplicity,  we  define 


(A. 3) 


is  tlie  Heaviside 
the  plasma- 


6♦^r)  =  Ef"^(r)  . 

Z  9  X 


(A.  4) 


Inside  the  electron  beam  (0  <  r  <  Rq)>  Eq.  (A.  3)  can  be  expressed  as' 


12 
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(A. 5) 


where 


.^2  -  2 
T  =  q 


2  /  2 
“pb^^b 


3_,2 


"  |[iii  -  (k  +  nkjj)V^  ■*■  l|k  ■^  nkjjlA/Yjlm] 


1>.  (A.6) 


Outside  the  electron  beam  (R.  <  r  <R  ),  Eq-  (A. 3)  reduces  to  the  free- 

U  C 

space  eigenvalue  equation 


The  solution  to  Eq.  (A.S)  that  remains  finite  at  r  ~  0  is 


6((.^^(r)  -  ijJ/TR)  .  0  r  <  Rjj  ,  (A. 8) 

where  Jj(x)  is  the  Bessel  function  of  the  first  kind  of  order  £,  and 
-  2  ? 

is  a  constant.  Noting  a  R«  >>  1  in  Eq.  (A.  2),  a'C  can  express  the 
^  n  U 

solution  to  Eq.  (A. 7)  as 


^♦^Ut<")  -  C(I7q„r)K7q„R,)  -  K7q„r)17q„R^)).  Rj,  <  r  <  R^,(A.9) 

where  and  are  modified  Bessel  functions  of  order  t  ,  and  C  is  a 
constant. 

The  boundary  conditions  on  6^^r)  at  the  surface  of  the  electron 
beam  are  given  by 

and 

from  Eq.  (A. 3).  Substituting  Eqs.  (A. 8)  and  (A. 9)  into  Eqs.  (A. 10)  and 
(A. 11)  gives 


A-2 
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TK 


0 


l.(q„) 


(A. 12) 

-  '’n-O  ’ 

where  the  "p*'i''ie"  notation  denotes  derivative  with  respect  to  the  complete 
argument  of  the  Bessel  function,  e.g.,  J’(TR,)  »  [dj  (x)/dx]  .  The 

£  \J  £  X*TKq 

expression  for  the  longitudinal  wave  admittance  h(q^)  in  Eq.  (A. 12) 

can  be  simplified  in  several  limiting  cases,  including  short  wavelength 
2  2 

perturbations  with  >>  this  case,  approximated  by 


h(qn)  =  -q„Rj;Cothq^(R^  -  R^)  ,  (A.  13) 

and  Eq.  (A. 12)  reduces  to 
J  * {TR  } 

J^(TRq)  '  ‘*n'^o‘^°'^‘’‘’i/’^c  ■  ^0^  •  (A.  14) 

Evidently,  the  right-hand  side  of  Eq.(A.14)is  a  very  large  positive 
number,  and  the  lovest-order  longitudinal  dispersion  relation  (for  A  -*■  0) 
can  be  approximated  by 


J^CTRq)  --  0  , 


(A. 15) 


where  T  is  defined  in  Eq.  (A. 6).  It  follows  from  Eq .  (A. 15)  that 


2  2  2 

T  Rq  »  6^  g.  .  s'  =  1,  2 . 


(A.  16) 


where  B  ,  is  the  s'th  aero  of  J  (x)  =  0.  In  this  regard,  Eqs.  (A. 8) 

£  •  S  £ 

and  (A. 9)  can  be  approximated  by 


64i^(r) 


otherwise  , 


(A. 17) 


where  ip  ,  is  a  constant.  Substituting  Eq.  (A. 6)  into  Eq.  (A. 16) 

£»8 


A3 
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and  making  use  of  Eq.  (A. 2),  we  obtain  the  longitudinal  dispersion  relation, 


2  2 


u  -  (k  +  nk„)V.  -H  i 
U  b 


|k  -f  nkn  I  a]  to  , 

- -  -Et  .  0 

3  2  • 

T'h 


(A.18) 


where  the  term  proportional  to  6^  has  been  neglected. 
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